Charge order induced by electron-lattice interaction in NaV 2 G*5 
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We present Density Matrix Renormalization Group calculations of the ground-state properties 
of quarter-filled ladders including static electron-lattice coupling. Isolated ladders and two coupled 
ladders are considered, with model parameters obtained from band-structure calculations for a'- 
NaV205. The relevant Holstein coupling to the lattice causes static out-of-plane lattice distortions, 
which appear concurrently with a charge-ordered state and which exhibit the same zigzag pattern 
observed in experiments. The inclusion of electron-lattice coupling drastically reduces the critical 
nearest-neighbor Coulomb repulsion V c needed to obtain the charge-ordered state. No spin gap is 
present in the ordered phase. The charge ordering is driven by the Coulomb repulsion and the 
electron-lattice interaction. With electron-lattice interaction, coupling two ladders has virtually no 
effect on V c or on the characteristics of the charge-ordered phase. At V = 0.46 eV, a value consistent 
with previous estimates, the lattice distortion, charge gap, charge order parameter, and the effective 
spin coupling are in good agreement with experimental data for NaV2 0s. 

PACS numbers: 71.10.Fd, 71.38.-k, 63.22.+m 



I. INTRODUCTION 

Since the observation of a phase transition at Too ~ 34 
K via magnetic susceptibility measurements/ the phys- 
ical properties of the low-dimensional inorganic com- 
pound NaV2 0g£ have been investigated intensively. The 
lattice distortions; 3 " the opening of a spin gap^ at Tqo 
or slightly below and the static charge disproportion 
6 between the V ions are the most remarkable proper- 
ties of this phase transition, which was at first identi- 
fied as a spin-Peierls transition like that in CuGe03. 7 
Later studies— found that zigzag charge ordering oc- 
curs below the phase transition in NaV20s, which is 
made up of quarter-filled ladders. It now appears that 
the opening of the spin gap may be induced primar- 
ily by charge transfes&Siifiiii in the low temperature 
superlatticeii2iiiiiii5ii& The discussion on the main driv- 
ing force of the transition is still going on. It was ar- 
gued recently that the Coulomb repulsion on an isolated 
ladder might be too small to cause charge ordering>iiii£ 
The inter-ladder coupling, which appears to cause the 
spin gap, would then also have to be responsible for 
the charge ordering. However, the electron-lattice cou- 
pling has been found to contribute to the charge-ordering 
transitioni 18 i 19 i 2n i 21 i 22 i 23 

In this work, we concentrate on the effects of electron- 
lattice interactions on the ground state of NaA^Os. 
Due to the asymmetric crystal environment of the V 
ions in NaV 2 5 , the strongest lattice coupling is a 
Holstein-like electron-phonon interaction involving the 
d xy electrons —2*24 We investigate the extended Hubbard 
model with Holstein coupling, introduced in Ref. l23l 
using high precision Density Matrix Renormalization 
Group (DMRG) calculations on isolated and on two 
coupled ladders. We show that, indeed, including the 
electron-lattice coupling causes charge ordering to appear 
on isolated as well as on two coupled ladders at signifi- 
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FIG. 1: Schematic depiction of the vanadium ladders 
in NaV2 05. The three types of hopping matrix ele- 
ments {t a ,tb,t a \,) and nearest-neighbor Coulomb repulsion 
(V a , Vb, Vab) are indicated. Two ladders are shown. 



cantly smaller Coulomb repulsions, with values which are 
consistent with most independent estimates. 

In Sec. |n] we present the model and discuss the pa- 
rameter values and observables. We discuss results for 
an isolated ladder in Sec. lIIII as a function of the nearest- 
neighbor Coulomb repulsion. Section Sec. II VI treats cou- 
pled ladders, and is followed by our conclusions. 



II. MODEL 

NaV205 is the only known quarter-filled ladder 
compound. 25 Ladder-like structures are formed by the 
vanadium ions and are only weakly coupled. A 
useful microscopic description is provided by an ex- 
tended Hubbard model (EHM) with Coulomb repul- 
sion between nearest-neighbor sites. The properties 
of the EHM without electron-lattice coupling were 
studied using mean field approaches ^SLSLS^ perturba- 
tion theory^ the Dynamical Mean-Field Theory^ ex- 
act diagonalization)2 3 j22i 3 i 3 ^ 3 .bosonization, 34 Quantum 
Monte Carlo^ and Cluster Perturbation Theory^ De- 
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tailed studies of the EHM were also performed using the 
DMRG)2LiS but not at specific values of the coupling 
parameters appropriate for NaV2C>5. 

A complete microscopic description of NaV20s must 
incorporate the lattice distortion observed by x-ray 
diffraction in the low temperature phased The distor- 
tions of the vanadium ions considerably change the dis- 
tance to the neighboring oxygen atoms directly above or 
below the vanadium site; 3 Recent calculations using Den- 
sity Functional Theory in the Local Density Approxima- 
tion (LDA) have shown that the main effect of these dis- 
tortions is a change of the local potentials, i.e., a Holstein- 
like coupling H e -im^ Indeed, the corresponding phonons 
are the strongest-coupling A g modest in NaVaOs; in 
the current study, we restrict the electron-lattice cou- 
pling to these distortions. The contribution of the lattice 
deformation to the Hamiltonian can be approximated by 
a parabolic potential Hi. This yields the microscopic 
Hamiltonian proposed and studied on small systems in 
Ref . El 

H = Hehm + Hi + H e -i (1) 

with 

Hehm = - E (4* c j<r + h - c -) 

+ U n^riii + VijUjUj, (2a) 
' («> 

Hi=k^2^, (2b) 

i 

H e -i = - C^ZiUi , (2c) 

i 

where n, = + n,-f is the occupation number and Zi 
the distortion at site i (in units of 0.05 A). The Hamil- 
tonian contains the effective lattice force constant k and 
a large Holstein constant C. The placement of the hop- 
ping terms t^j and the Coulomb repulsion in the lat- 
tice structure are depicted in Fig. ^ The parameters are 
taken from LDA calculations 24 ; the hopping amplitudes 
tij are compatible with earlier calculationsiS^S The lat- 
tice parameters C and k were extracted by comparing 
the total energy and the inter-ionic forces in distorted 
and undistorted lattices. They correspond to the V-O- 
stretching perpendicular to the ladder plane observed in 
the 970 cm -1 phonon mode. As in Ref. |23L we take 

t a = 0.35 cV = 2 t b , t ab = 0.17 t a and 0.33 t a , 
V a = V b = V ab = V, U = 8.0t a , 

k = 0.125 t a , and C = 0.35t a . (3) 

We compare results for two different values, 0.17£ a and 
0.33 t a , of the interladder hopping because existing es- 
timates of these parameters differ i 2 ^ 2 ^* 3 ^ The Coulomb 
repulsion is very difficult to compute; estimates in the lit- 
erature vary stronglyAiLSLSiSMSi In our study, we take 
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FIG. 2: Square of the charge order parameter, m,Q , on an 
isolated ladder without coupling to the lattice, as a function of 
the nearest neighbor Coulomb repulsion V for several ladder 
lengths L. The inset shows the finite size scaling data collapse, 
with exponents v = 1, /3 = g, and critical Coulomb repulsion 
Vc = 2.31(2). 

V a = V = V a b = V for simplicity, and investigate the 
model as a function of V. 

We have performed DMRG calculations on isolated 
ladders of length of up to L — 80 and on two coupled 
ladders of length of up to L — 24, applying open bound- 
ary conditions along the chains (6-direction) . In the case 
of two coupled ladders, we have taken periodic boundary 
conditions in the a direction. 

We have measured the charge order parameter mco> 
defined by 

ij 

where Q = (jr, 7r) and N is the total number of sites in the 
system. The occupation number of the two inequivalcnt 
V ions is then (1 ± mco)/2. 

The charge gap Ac(L) and the spin gap Ag(L) are 
determined using 

A C (L) = ^[E a (L,N + 2) + E (L,N-2)-2E a (L 1 N)} 

A S (L) = E (L, N, S z = 1) - E (L, N, S 2 = 0), (5) 

where Eq(L, N) is the ground state energy of the system 
with L rungs per ladder and N electrons^ Results are 
extrapolated to L — oo using linear and quadratic fits in 
1/L and by applying finite-size scaling. 

III. ISOLATED LADDER 

In this section, we present DMRG calculations on iso- 
lated quarter-filled ladders (V a b — 0, t a b = 0). We con- 
sider systems with (C = 0.35) and without (C = 0) 
electron-lattice interactions and compare the results. 

Isolated quarter-filled ladders without coupling to the 
lattice were investigated by Vojta et. alM^ for a range 
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FIG. 3: Optimal zigzag distortion z op t, extrapolated to L — 
oo, as a function of V on an isolated ladder system with 
electron-lattice coupling. The inset shows the linear 1/L ex- 
trapolation for some values of V. The experimentally de- 
termined distortion in NaV20s corresponds to V* = 1.3 t a , 
indicated by the dotted lines. 



of parameters. The parameters studied did not, how- 
ever, include the couplings relevant for NaV20s, which is 
characterized by large hopping along the rungs, t a ~ 2tb, 
large Coulomb repulsion U ~ 8t Q , and charge ordering. 
Ref. IH found qualitatively different behavior for t a < tb, 
where the spin gap was found to be finite, and t a > tb, 
where the spin gap was found to vanish. 

In Fig. [21 we display the square of the charge order 
parameter Wq as a function of the nearest-neighbor 
Coulomb repulsion V for different system sizes. At a 
critical value V c above 2t a , there is a quantum phase 
transition to a phase with finite charge order. The essen- 
tial features of this transition can be described by a model 
with a single charge degree of freedom on each rung, i.e., 
a pseudospin. When the hopping between rungs is ne- 
glected, one arrives at the Ising model in a transverse 
field (IMTF)^ which can be solved exactlyii Indeed, 
a finite-size scaling analysis for m co in which L 2 ^' u m co 
is plotted as a function of L x l v Qy — 1) with the critical 
exponents v = 1, (3 = | of the IMTF, shown in the in- 
set, collapses all data points onto a universal curve and 
yields a critical Coulomb repulsion of Vc = 2.31(2) t a . 
We have also performed a similar scaling analysis with- 
out fixing v and (3, and obtain a reasonable data collapse 
with v = 1.0(1), [3 = 0.125(20), and V c = 2.31(4). We 
note that (3 = | (which is also the exponent of the 2d 
Ising model) is fairly close to the value extracted from 
experimental results for NaV205f£^, about 0.15-0.20. 
However, the effective Coulomb coupling V in NaV2 0s 
is probably considerably below V c ^ 2.31 £ Ql 2JL28i22^SUtt 
Therefore, the Coulomb repulsion alone is insufficient to 
cause the charge ordering in NaX-^Os. 

When the electron-lattice coupling is taken into ac- 
count, the behavior of the charge order parameter 
changes drastically, as we will now show. The ground 
state energy E ({zi}) of the Hamiltonian H (Eq. QJ) i s 
a function of the independent classical lattice distortions 
Zi on each site. First, we have determined the optimal 
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FIG. 4: Square of the charge order parameter, m-co, calcu- 
lated at the optimal distortion z op t on an isolated ladder with 
electron-lattice coupling as a function of V and for several 
ladder lengths L. The solid line is the result of a linear 1/L 
extrapolation, illustrated in the inset for selected values of V. 
The dotted lines mark V* = 1.2>t a - 

distortion pattern on ladders of up to 16 sites using clas- 
sical Monte Carlo simulations carried out in the space of 
all Zi. 

In these simulations, the ground state energies are 
determined by exact diagonalization. We use parallel 
tempering 4 - to find all relevant distortion patterns. We 
find that at low V the optimum lattice distortion (lowest 
total energy) is Zj = 0, i.e., no distortion. Above a criti- 
cal value V c w 0.95 t a , there are two degenerate optimal 
configurations, with finite zigzag lattice distortions 

z i = ze i * H * ) Q=(7r,7r), (6) 

where Ri labels the lattice sites. The zigzag pattern of 
the lattice distortions agrees with the experimentally ob- 
served pattern^ In the following, we therefore assume a 
zigzag pattern. We determine the ground state by min- 
imizing Eq(z) as a function of the lattice distortion z, 
similarly to Ref. I23I The position of the minimum de- 
fines the optimal distortion z opt . 

The optimal distortion z opt and the extrapolation to 
the thermodynamic limit are presented in Fig. The 
distortion becomes finite at a critical Coulomb repulsion 
Vc = 0.95(1). From Fig. [3] and from the experimentally 
determined size of the zigzag-distortion^ (z exp rs 0.85 x 
0.5 A), we obtain an estimate for the effective Coulomb 
repulsion in NaX^Os of V* = 1.3 t a , well within the range 
of most earlier estimatea 17 i 28 i 29 i 30 i 40 . 

The square of the charge order parameter m^, (L), 
calculated at the optimal distortions z op t{L), is shown 
in Fig. along with a finite-size extrapolation. Order 
sets in at the same V as the lattice distortion z opt . For 
the 1/L-extrapolated values, m co is proportional to the 
square of the optimal distortion z^ pt at all V. In a finite- 
size scaling analysis (not shown), the scaling of mco is 
consistent with the mean- field exponent (3 = ^, but not 
with the IMTF exponent j3 = |. 

A comparison of Fig. |3 and Fig. Q] illustrates the sub- 
stantial decrease of the critical Coulomb repulsion due 
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5: Charge gap Ac on an isolated ladder as a function 
with (solid) and without (dashed) electron-lattice cou- 
pling. The results at V* = 1.3 t a are indicated by dotted lines. 
Inset: linear finite-size extrapolation in 1/L for systems with 
(dist.) and without (no dist.) lattice distortion. 
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FIG. 7: Effective magnetic exchange coupling J e ff (Eq. Q) 
as a function of inverse lattice size 1/L, for different values of 
V, with and without lattice distortion. 



to the electron-lattice coupling. The critical value V c de- 
creases from V c = 2.31(2) t a without electron-lattice cou- 
pling to V c — 0.95(1) t a with electron-lattice coupling. 
The critical exponent of mco is clearly smaller than 
unity in both cases, but changes from /3 ~ g to (3 w | 
when the electron-lattice coupling is switched on. At the 
coupling V* at which the lattice distortion matches the 
experiment value, we find m^. ss 0.37, very close to the 
experimental results of about 0.35i£ and 0.37^ 

The charge gap Ac is also influenced by the presence of 
electron-lattice coupling, as shown in Fig. |S| The charge 
gap increases with increasing V and lattice distortion and 
does not vanish^ at V — 0. At V* — 1.3 t a , the charge 
gap agrees reasonably well with the experimental value 
for the optical gap in NaV205^ 

The spin gap of an isolated quarter-filled ladder with- 
out electron-lattice coupling has been shown by Vojta et 
al~ to vanish at moderate V when the rung hopping t a 
is larger than tb, as is the case in NaV2C>5. In Fig.^a), 
we show some representative results of our calculations 
of the spin gap, with and without electron-lattice cou- 
pling, together with a quadratic 1/L-extrapolation. We 
find that the spin gap extrapolates to zero on an isolated 



ladder at all V in both cases. Charge ordering on an iso- 
lated ladder is therefore not sufficient to induce a spin 
gap. 

The magnetic behavior of the quarter-filled Hubbard 
ladder can be approximated by an antifcrromagnetic 
Heisenberg modeli^^iiSiiS We determine an effective mag- 
netic exchange interaction J c g by equating the finite-size 
spin gap in our model to that of a Heisenberg chain with 
exchange constant J, i.e., 

A^ ciscnbcrg (L) A" ubbard (L) 

J = J oS (L) ■ (7) 

The results are plotted in Fig. as a function of inverse 
system size 1/L. There is only a weak L-dependence, 
and the scaling is linear in 1/L. Previous work using 
exact diagonalization 23 showed that the behavior of the 
dynamical spin correlations S(k,uj) for the EHM with 
electron-lattice coupling also closely resemble those of the 
Heisenberg chain. Remarkably, our results for J g agree 
very well with the estimates from the spin dispersion in 
Ref.lH 

The resulting values of J e g are plotted in Fig. [S] as 
a function of V. As F is increased, the charge order 
parameter increases, and J g decreases. This behav- 
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FIG. 8: Effective magnetic exchange coupling J c s on an iso- 

FIG. 6: Spin gap A s for (a) an isolated ladder and (b) two lated ladder as a function of V, with (solid line) and without 

coupled ladders. Lines represent a quadratic fit in 1/L. The (dashed line) electron-lattice coupling. The dashed-dotted 

boxed descriptions of the curves are in the same vertical order line is the result of perturbation theory 38 for the undistorted 

as the data. lattice, valid only at large charge order. 
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FIG. 9: Square of the charge order parameter, rn% , for two 
coupled ladders (4 x L, t a t = 0.33 t a ) with periodic boundary 
conditions in the a direction and with and without lattice 
distortions. The results for infinite length are obtained using a 
linear extrapolation in 1/L. Also included for comparison are 
similarly extrapolated results for an isolated ladder (2 x oo). 

ior is in accordance with the experimental observation 
that Jeff becomes smaller at lower temperature, where 
charge order increases. Part of the reduction in J c ff is 
due directly to the charge occupation of neighboring sites 
along the chains, which implies a reduction by a factor 
of (1 — m^ )£ Our results, however, are far from such a 
simple quadratic dependence on rr? co . 

IV. COUPLED LADDERS 

The occurrence of a spin gap in NaV2C>5 may largely be 
due to the coupling of laddersi 8 i 9 i 10 i n We have therefore 
studied a system of two coupled ladders with periodic 
boundary condition in the a direction. We have deter- 
mined the pattern of optimal lattice distortions in the 
same way as for the isolated ladder (for t ao = 0.33 t a ), 
using exact diagonalization on a 16-site system, and find 
a simultaneous zigzag distortion on both ladders to be 
optimal. Due to the lattice structure (see Fig. QJ, the 
energy is invariant under a shift of the distortion by one 
rung on either ladder. Using the DMRG, we have then 
examined larger systems. We have determined the opti- 
mal distortion z op t from the minimum of the ground state 
energy Eq(z) as a function of the uniform zigzag distor- 
tion. The result is very similar to the lattice distortion in 
the isolated ladder at the same V. The effective Coulomb 
repulsion determined from the experimental lattice dis- 
tortion is therefore still V* = 1.3 t a . 

Without electron-lattice coupling, convergence of the 
DMRG calculations is difficult to achieve, restricting the 
largest length to at most L = 12 at large V. The three 
right-most curves in Fig.^Jshow the results for the charge 
order parameter as a function of V for lattice sizes L = 8, 
L = 12, and for a 1/L extrapolation. For comparison, 
the results for a single ladder without electron-lattice 
coupling extrapolated in the same way are also shown 
(middle solid line marked with stars). Clearly, the sys- 
tem consisting of two coupled ladders tends to order at 



even larger values of V than the isolated ladder, leading 
to an even larger discrepancy with the estimates for V in 
NaV 2 5 . 

At the lattice sizes we were able to reach, the effects 
of the open boundaries, especially effects caused by the 
different number of neighbors at the boundary, are con- 
siderable. While varying the boundaries, e.g., compen- 
sating for missing neighbors, did affect the results shown 
in Fig. 1^1 it did not change the overall tendency to shift 
the phase transition to larger V. 

For systems with electron-lattice coupling, we were 
able to reach larger sizes, up to 4 x 16 sites at large V 
and up to 4 x 24 at V* — 1.3 t a . The results, as can be 
seen in Fig. [51 are almost identical to the results on a 
single ladder of the same length L, and are almost inde- 
pendent of length in the ordered phase. Charge ordering 
still takes place at V ~ 0.95 t a . 

The spin gap of the coupled ladder system is displayed 
in Fig.HJb) for V* — 1.3 t a . Since estimates of the inter- 
ladder hopping t ao in NaA^Os varyj2i2iiiSi we have per- 
formed our calculations for two values, t ao — 0.17 t a and 
0.33 t a . At t a b — 0.17 t a , the spin gap for two coupled lad- 
ders both with and without lattice distortion is similar in 
value and finite-size behavior to that for the isolated lad- 
der, which vanishes in the thermodynamic limit. When 
the inter-ladder hopping is increased to t ao = 0.33 t a , the 
spin gap becomes smaller and its finite-size dependence 
appears to change in both cases. The physics of the indi- 
vidual ladders thus appears to be noticeably affected by 
inter-ladder coupling at this larger hopping. 

Extrapolations of the gaps using a fit to a quadratic 
polynomial in 1/L arc indicated by the lines in the fig- 
ure. The value of the spin gap in the L — » oo limit 
for all parameter values is approximately zero, |Ag(_L = 
oo)| < 0.01 t a . Some of the extrapolations in Fig. 
yield values that are slightly negative. The magnitude 
of these negative values gives as a minimal error esti- 
mate for the extrapolation. Since the magnitudes of the 
positive extrapolations are of similar size, we conclude 
that the extrapolated spin gaps are zero to within er- 
rors. In addition, the small extrapolated gaps are signif- 
icantly smaller than the spin gap found experimentally 
in NaV 2 5 , A s » lOmeV w 0.03 t a . The spin gaps at 
other values ofV (V — 0.9 t a , 1.5 t a ) behave similarly. 

We conclude that when two ladders are coupled, there 
is no significant indication that a spin gap opens. This re- 
sult is not unexpected, however, because the mechanism 
which has been proposed as the cause for the spin gap 
in the charge-ordered systemj&Siifiiii namely an in- ladder 
dimerization of effective spin couplings, requires a charge- 
ordering pattern with a period of four ladders in the a 
direction; such a pattern is observed in NaV205ii2ii&±2ii£ 



V. CONCLUSIONS 

In this paper, we have investigated the influence of 
Coulomb repulsion and electron-lattice interactions in 
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EHM+lattice 


Experimental 




V = 1.3 t n 


data 


Distortion/0.05A 


0.84 


0.85, Ref.£ 


m 2 co 


0.37 


0.35 - 0.37, Refs. 10.14 


Jeff 


63meV 


60meV, Ref.^0 


Ac 


0.7 eV 


0.9 eV, Ref. 44 



TABLE I: Comparison of experimental data for NaV^Os to 
the results for the Extended Hubbard Model with lattice cou- 
pling, at V* = 1.3 t a - 



quarter-filled ladder materials such as NaA^Os by treat- 
ing the extended Hubbard model with and without cou- 
pling to static lattice distortions in the c direction. 
Our calculations show that the electron-lattice coupling 
drastically affects the physical properties of the low- 
temperature phase. It causes lattice distortions to appear 
concurrently with charge ordering and with the same spa- 
tial zigzag pattern. The transition to the charge-ordered 
phase is shifted to a much smaller critical value of the 
nearest-neighbor Coulomb repulsion (V c ~ 0.95 t a ). The 
charge gap is increased and the effective magnetic ex- 
change is reduced by the electron-lattice coupling. These 
results remain unchanged when two ladders are coupled. 
The spin gap extrapolates to zero in all cases. 

When the electron-lattice coupling is included, the 
properties of both the isolated ladder and two coupled 
ladders are in remarkably good agreement with exper- 
imental data for NaV20s, at an effective Coulomb re- 



pulsion V* — 1.3 t a determined by matching the experi- 
mental lattice distortion. This Coulomb repulsion is well 
within the range previously estimated for NaV205. As 
shown in Tabled we find simultaneous agreement for the 
zigzag distortion, the size of the charge order parameter, 
the effective spin interaction, and for the charge gap with 
experimental data. 

While results for the isolated ladder without electron- 
lattice coupling would also be reasonably consistent with 
the experimental values, an unrealistically large value for 
the Coulomb repulsion, V ~ 2.5 t a , would be required, 
and the experimentally observed lattice distortion would 
be not present. The value of the Coulomb repulsion nec- 
essary for such an agreement becomes still larger when 
two ladders without electron-lattice interaction are cou- 
pled. 

We conclude that an interplay of charge ordering 
and lattice distortion can drive the phase transition 
to a charge-ordered state in NaV205, independent of 
the occurrence of a spin gap. The explanation of the 
spin gap will likely require a unit cell with at least four 
coupled ladders. Corresponding work is in progress. 
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